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Abstract. We study the discrete version of the p-Laplacian. Based on its variational properties we discuss some 
features of the associated parabolic problem. Our approach allows us in turn to obtain interesting information 
about positivity and comparison principles as well as compatibility with the symmetries of the graph. We 
conclude briefly discussing the variational properties of a handful of nonlinear generalized Laplacians appearing 
in different parabolic equations. 

1. Introduction 

The discrete Laplacian is a well-known object in graph theory. In his seminal investigation on electric cir- 
cuits [2^ , Kirchhoff introduced it as 

where D is the (diagonal) degree matrix and A the adjacency matrix of an undirected graph G. He then went 
on to observe that for any orientation of the graph the associated (signed) incidence matrix I satisfies 

(1.1) A=IZ^, 

that this representation allows us to prove that all the minors of A have the same value, and that this value is a 
natural number that in fact agrees with the number of spanning trees of G; moreover, is always an eigenvalue 
whose multiplicity coincides with the number of connected components of G. 

A possible reason for considering A as a discrete analogous of the usual Laplace operator is that, if one sets 
up a system of linear equations for the potential / at each node v of an electric circuit in accordance with the 
(linear) laws of Kirchhoff and Ohm, then 

A/ = 0, 

i.e., / satisfies a Laplace- type equation like its continuous pendant; this is well explained in [lil Chapt. 1]. 
Also after Kirchhoff's investigation, A has proved a remarkable object in algebraic graph theory. The role of 
the second smallest eigenvalue of A in relation with the study of connectivity properties of a graph has been 
emphasized since the 1970s by Fiedler and others, cf. the survey article [s^l- Moreover, it is nowadays known (see 
e.g. [131) t^^t discrete Laplacian (possibly after a suitable re-normalization) is tightly related to the Dirichlet- 
to- Neumann (i.e., the voltage-to-current) operator of the metric structure associated with the discrete graph - 
the so-called quantum graph. A thorough theory of linear discrete elliptic operators is nowadays available. This 
field has experienced successful interactions between the communities of researchers working on graph theory, 
functional analysis and potential theory, see e.g. [27l.[30l|. 

In the 1960s and 1970s, pioneering investigations by (among others) Minty, Rockafellar, and Zemanian aroused 
broad interest in the theory of nonlinear electric circuits. In those years also the theory of monotone operators and 
subdifferentials was being substantially developed. It is therefore no surprise that many investigations focused 
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on those nonlinear electric circuits whose associated voltage-to-current operator is maximal monotone. In the 
same years, the (continuous) p-Laplace operator Ap, i.e., the subdifferential of the energy functional 



(1.2) 



LP{n) ^||Vw||^p G [0,oo] 



began to be studied by Aronsson, DiBencdctto, Ural'ccva, and many others (cf. llTll 
survey on this topic). It was in this cultural climate that Yamasaki proposed in [44 
integral of order p G (1, oo) on a graph G with node set V: by this the mappingQ 



for an introduction and 
1 the notion of Dirichlet 



Sp : R 



V 



P 



E 

v,wGV 



|/(v)-/(w)r e [0,oo], 



is meant. It seems that Yamasaki's definition went forgotten until the early 1990s, when (op reappeared as part 
of the final remark of a short note (39| on linear potential theory by Soardi. Ever since, many authors have 
discussed properties of this energy functional and of its subdifferential, the discrete p-Laplacian defined by 



Ap/(v) := ^ |/(v) - /(w)r-2(/(v) - /(w)). 



V e V. 



(Observe that for p = 2 we recover the discrete Laplacian; however the discrete Laplacian is positive definite 
whereas the common Laplace operator is negative definite). 

Most investigations have been devoted to potential theory and eigenvalue problems: we mention, among 
others, P, [lH, [s^l and a long series of articles by Agarwal, O'Regan and coauthors that begins with [2^. In 
recent years, several works have emphasized the role of the discrete p-Laplacian in spectral clustering and image 
processing, cf. 



41|,|45|. 

The introduction of (op is motivated by the fact that the finite difference 

Z^/(v,w) :=/(v)-/(w) 

can be seen as a discrete analogue of (minus) a directional derivative - this fact is classically used in numerical 
analysis, giving rise to the notion of backward difference operator. If we impose an orientation on a graph G and 
regard functions from V to R as difference 0-forms (or 0-chains) and functions from E to R as difference 1-forms 



(or 1-chains), the corresponding derivation is exactly 
incidence matrix 



— J> 



", the transpose of the oriented |V 



(1.3) 

where 



(t+) and X := (tve) ^'^^ defined by 

1 if V is initial endpoint of e, 
otherwise. 



if V is terminal endpoint of e, 
otherwise, 



respectively. Hence, Yamasaki's functional S'p can be written as 

' 11^^/11?. 



S'r, 



3/^ 



V 



(E) 



G [0,00], 



and 

(1.4) : 3 / ^ X{\X^ fr^X^ f) G R^, 

in apparent analogy with the continuous p-Laplacian operator and its associated energy functional. 
Just like in the continuous case, where (after applying the divergence theorem) the conservation law 

d 
It 



(j){t, x)dx 



diMi{t,x)dx / f{t,x)dx 



^ Throughout this paper wc write v ~ w if there is an edge between v and w. 
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is considered, in the discrete case we have the conservation law 

I J2 Ht, v) = - Y.(^j)it, e)+J2 fit, v) 
vsv vev vev 

and Pick's law 

j-(i,e) = -c(I^0)(t,e) 

for a general function c can be used to derive a general (time-continuous, space-discrete) differential equation 
governing a flow on G, cf. 0, § 2.5.5]. Indeed, for c = 1 we obtain exactly the linear heat equation for the 
discrete Laplacian. Considering different terms c we obtain different equations. In particular, opting for the 
nonlinear Darcy's law 

c = c(t,e,(I^0)(i,e)) |X^<^(t, e)^-^, 

we end up with 

I ^ cf,{t, v) = -Y,m^<i>r'i^m. v) + ^ fit, v, ^(t, v)). 

vsv vev vsv 

In view of ()I.4p . this is the summed version of the parabolic problem associated with the discrete p-Laplacian 
Ap - let us call it the discrete p-heat equation. Observe that, unlike its continuous counterpart, it is a backward 
evolution equation. 

The present paper is devoted to study the basic features of this equation and is structured as follows. After an 
introductory section where the relevant functional setting is introduced, we proceed by formulating in Section [3] 
our main well-posedness result for the discrete p-heat equation. A Galerkin-type method can be applied to 
yield a solution and, unlike in the usual context of PDEs on domains, we can even interpret the converging 
Galerkin sequence yielded as a sequence of solutions to the same equation on induced subgraphs (i.e., on growing 
subsystems). Using the theory of nonlinear Dirichlet forms, we can show that these solutions arc associated with 
semigroup of nonlinear contraction on all ^''-spaces, 1 < <? < oo, as well as on cq. This is done by discussing 
the invariance property of several relevant closed convex subsets under the semigroup generated by —r^p. We 
show that the nonlinear semigroup generated by the discrete p-Laplacian consists of irreducible, sub-Markovian 
operators. In the special case of p = 2, these results had already been obtained in [23|. 

Similar methods are applied in Section 0] to prove how these semigroups interact with the symmetries of the 
graph. This is more delicate than in the continuous case, since the discrete p-Laplacian is not a local operator. 

Finally, in the Appendix we overview a few popular generalizations of the discrete Laplacian, propose some 
more and suggest how the variational structure of all can be used to discuss a broad class of discrete nonlinear 
diffusion-type problems, including a discretized porous medium equation. 

Because the discrete p-heat equation is simply a dynamical system, most of our results follow from the standard 
theory of ordinary differential ecruations whenever G is a finite graph. The actually interesting case is hence that 
of infinite graphs. We refer to [15|, Chapter 8] for an introduction to infinite graph theory. On the analytical 
side, standard references for nonlinear Cauchy problems on Banach spaces are the monographs [1, [s^ . 

2. General setting and the energy functional 

We consider throughout a (finite or countable) set V and a (non-symmetric) relation 

E C { (v, w) e V X V : V ^ w} 

such that for any two elements v, w G V at most one of the pairs (v, w), (w, v) belongs to E. We refer to the 
elements of V and E as nodes and edges, respectively, and regard G as a graph whose incidence matrix is I 
introduced in (jl.3p . If e = (v,w), we denote the initial and terminal endpoint of e by 

eo := V and ei w, respectively, where e ~ (v,w). 

Inspired by [ivl § 1], throughout this article we consider G as a weighted graph, i.e., as a quadruple (V, E, /i, v) 
where /z : E — >• R+ and : V — )• IR+ (the unweighted case corresponds to /it = 1 and v = 1). 
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Remark 2.1. Actually, by construction G is oriented. This is reflected already in the definition of the (signed) 
incidence matrix I. However, orientation is only a technical tool that will make some parametrizations easier: 
Unless otherwise underlined, throughout most of this paper only the absolute value of the numbers 

vev 

will play a role. Thus, all our results ( with the exception of those in Section are still valid if any of the edges 
is given the opposite orientation. 

We have already introduced the (possibly infinite) incidence matrix 

X : ^ and its transposed : ^ R^. 

Along with them, the structure of the undirected graph associated with G is mirrored by the |V| x |V| adjacency 
matrix A :— (ovw) defined by 



(2.1) ttv 



/.t(e) if e = (v,w) or e — {w,y) e E, 
otherwise. 



Definition 2.2. A graph G := (V, E, fi, v) is called outward locally finite if 

deg^(v) := ^ tJeM(6) < for allv e\/ and some > 0. 

eeE 

It is called inward locally finite if 

dcg^(v) := /,^^/i(e) < M~ for a llv ^\/ and some > 0. 

eeE 

It is called locally finite if it is both inward and outward locally finite, i.e., if 

dcg(v) := dcg^(v) + dcg^(v) < My for all v £\/ and some My > 0. 

J/deg'*' £ 0{v), i.e., if there exists > s.t. 

dcg+(v) < A/+j/(v) for all v G V, 

then G is called outward uniformly locally finite. We define inward uniform local finitencss and uniform local 
finiteness likewise. 

We say that G has finite surface and finite volume if 

|V|^ := ||i^llfi(V) < oo and |E|^ := ||^||fi(E) < oo, 

respectively. 

Observe that a graph with finite volume is necessarily uniformly locally finite. 

Assumptions 2.3. Throughout this paper we impose the following assumptions on G := (V, E,/i, :/). 

• /i(e) > for all e € E. 

• v{y) > for all v e V. 

• G is locally finite. 



The weights n, v define in a natural way a metric and measure structure on G. As we will see in Lemma 13.51 
the (discrete) p-Laplace-Beltrami operator associated with these coefficients is given by 

(2.2) Ap/ := -X{pi\X^ f\v-^X^ f), f e . 

u 

However, we would like to allow for a general class of (non-degenerate) elliptic operators, possibly in non- 
divergence form. To this aim, we will have to consider general weights a and d. In order to avoid degeneracies, 
certain compatibility conditions should hold. 
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Assumptions 2.4. Throughout this paper the functions a : E — )- R and d : V — )- R are assumed to satisfy the 
following conditions. 

• There exist k, K > s.t. 

Ka{e) < fi{e) < Ka{e) for all e G E. 

• There exist 9,Q > s.t. 

ed{\/) < i/(v) < ed{\/) for all V e V. 

Clearly, if there exist (, Z > s.t. 

C < ti-< Z (resp., s.t. C < Z), 

then one can choose a = 1 (resp., d = 1). 

For p G [1, oo] we will consider the weighted sequence spaces £^{E) and £'^{V) defined by 



along with 

In other words, 



J2 HeWaie) and ||/||^,(^) := ^ |/(v)|''d(v), 

eGE vGV 



<"/(V) :={/G£?(V):I^/G£?;(E)}. 
If G £l{y) : ||T^/||Pp = ^a(e)|/(eo) - f{e,)\P < +oo 

I eGE 



Remark 2.5. Letp^ [l,oo]. Observe that another way of expressing uniform local finiteness of G , i.e., thatdeg G 
0{v), is saying that ^^gg(V) is continuously embedded in ^^(V). Also observe that under the Assumptions \2.4\ the 
weighted norms \\ ■ ||^p and \\ ■ jj^p, and also \\ ■ jj^p and \\ ■ jj^p are equivalent. In particular, wj^'P;'^ and w^'^''^ agree. 

In the theory of discrete calculus, the rule of thumb is to perform the following substitutions: 

• scalar functions — ?> vectors of the node space, 

• vector fields — >■ vectors of the edge space, 

• gradient of a scalar function at a point x — > evaluation of l'^(f> at an edge e, 

• divergence of a vector field tp at a, point x — >■ evaluation of lip at a node v. 

This interplay between continuous and discrete setting is evocative and already thoroughly investigated (see 
e.g. (l^l). Actually, '"^^'^'^(V) will play the role of a weighted Sobolev space. 

Lemma 2.6. For all p, q G [1, oo), £a{^) o,nd ^^(E) are separable Banach spaces, and so is u'^'2'^(V) with respect 
to the norm defined by 



For all p G [l,oo], u'^'d'^(V) is continuously and densely embedded into i'^(V). If moreover p G [l,oo), then 
^a'd'^(^) separable. If p G (l,oo), then w^'^'^(V) is uniformly convex (and hence reflexive). 

Proof. Separability of ^^(V) and follows from countability of V, hence of E. Considering 

<d'(v)3/^(/,i'^/)e^^(v)x^?;(E), 

which is an isometry if the Cartesian product on the right is endowed with the ^^-norm, shows that w^'^'^ is 
separable for all p G [1,cxd), and uniformly convex if p G (1,cxd). Furthermore, since the space coo(V) C R^ of 
functions with finite support is dense in ^^(V), so is w^'^' (V). □ 
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Lemma 2.7. Under our standing Assumptions fXM let p G [l,oo]. If G is outward uniformly (resp., inward 
uniformly, uniformly) locally finite, then (resp.,X~,X) is bounded from £p{E) to £^(V). The converse impli- 
cation holds for p G [1, oo) . It also holds for p ~ oo if additionally there exist /i, i> > s.t. /i(e) < // and v < 
for all e G E and v G V. 

In particular, G is uniformly locally finite if and only if ^^'^'^(V) = ^d(^) /''^ '^^^ P ^ [^i 

This result is a generalization of 0, Prop. 6] to weighted graphs. 
Proof. Take / : V — > IR and observe that 

E^'el/Wr = l/(eo)r - max4|/(v)|P for all e G E and ah p G [l,oo), 

vGV 

because there exists exactly one v G V s.t. /,+ 7^ 0. 

Take first p G [f , 00). By Fubini's theorem, one has for all / : E — IR 



eeE 

eeE VveV / 



vGV eeE 

El/(v)Reg+(v)=. 11/11^, . 



vev 

This shows that 1+^ is an isometry from ^^(E) to ^^cg(^)' hence by Remark 12.51 it is bounded from £f^{E) to 
£^(V) if and only if G is uniformly locally finite. One concludes that 1+^ is is bounded from ^^(E) to ^^(V) if 
and only if G is uniformly locally finite. 
For p — 00 one has 

= max|/(eo)|M(e) 

^ eGE 

= maxmaxtjg|/(v)|/i(e) 
eeE vev 

- E"^^?^^^el/MlM(e) 
^ — ^ vev 

eeE 

= max|/(v)|E'-JeA*(e) 
vev ^ — ^ 

eeE 

= max|/(v)|deg(v) = ||/||^~ 
vev ^ 

To prove the converse implication, take a sequence (v„)„giM C V s.t. 

nv{\/n) < deg+(v„) 

and consider the functions ?i„ : E ^> R defined for all 71 G IN by 

M„(e) := 1^+, 



PARABOLIC THEORY OF THE DISCRETE p-LAPLACE OPERATOR 



7 



where E+ ::={eG E;i+g7^0}, the set of edges outgoing from v„. Then ||m„||£oo < fi for ah n G N, but 



max 



eeE 



> J2 (^leA^(e)) ^ 

= deg+(v„)^ 

> n — , 

hence Um„_^oo H^^^Wn ||foo = +c)o. This concludes the proof. □ 

Remark 2.8. (1) Recall that the oriented incidence matrix I of a finite graph G and k connected components 
has rank |V| — k, see e.g. f2^ . § 8.3.1]. Hence, X is never surjective (resp., is never injective). On the other 
hand, I is injective (resp., 2^ is surjective) if and only if it has rank \ E\, which is the case exactly when G is a 
forest. It is easy to see that, more generally, 1^ is surjective whenever G is a uniformly locally finite forest. 

(2) Let now G he infinite and connected. If f G R^, f ^ 0, satisfies f = 0, then necessarily f is constant 
and hence does not belong to unless \\M\ep < oo, i.e., unless G has finite surface (for p < oo) or unless v 

is hounded (for p = oo). 

Taking into account Remark 12. 8[ a direct computation yields the following. 

Corollary 2.9. Let p G [1,cxd). Then a necessary condition for the operator 2^ ; ^^(V) —> ^^(E) to he an 
isomorphism is that G he infinite; a sufficient one is that G he a uniformly locally finite tree with infinite surface. 

Following 44 1, we consider an energy functional 

4, : / = i^a(e)|/(eo) - /(ei)|^ / G R\ 

P " 

eeE 

The definition of £'p is clearly independent of the orientation of G. While by allowing the value S'p{f) = +oo we 
can (and indeed shall) consider (^p as defined on the whole space ^^(V), one also considers the effective domain 

D{Sp) {/ G ^^(V) : Sp{f) < +^}, 

which is clearly non-empty: one says that (op is proper. In fact, D[S'p) = u'ali'^(V). In view of the mentioned 
analogies between discrete and continuous calculus ^p is the discrete companion to the functional introduced 
in (lL2l). 



Remark 2.10. A natural analogue of our setting arises in the theory of electric circuits, following the scheme 
• m(v) voltage at v, 



a(e) — >■ conductance of the hranch corresponding to e. 



• l^u{e) — > potential difference between the endpoints of e, 
• 

• a{e)X^u{e) current flowing between the endpoints of e, 

• setting Dirichlet boundary conditions on a certain set Vq of nodes — >■ grounding the nodes in Vq , 

• replacing the values of a function u in two nodes v, w G V by their average — >■ shorting the nodes of the 
corresponding electric circuit, 

• replacing the value of a in an edge e by ^ cutting the branch of the corresponding electric circuit. 



^ For the sake of notational simplicity, we stipulate that ifi^) will be written as (indifferently) either ^/'(v, w) or </'(w, v), whenever 
e is an edge with endpoints v, w and ip is a function from E to R. 
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Observe for example that cutting branches induces a shift in the Rayleigh-type quotient 

11/11^/ 

'■d 

which is the natural object to minimize if one looks for eigenvalues, cf. 0/ (i.e., for those A 6 IR s.t. 

A|</.(v)r2(^(v) ^/IpV'W, vgV, 

has a non-trivial solution ip). This shift is in accordance with the intuition that diffusion is slower on a graph 
with fewer edges and corresponds to Lord Rayleigh's Monotonicity Law (fldl - § 1-4])- will not go into details 
of this interplay and refer instead to to the beautiful booklet fldl l. where the linear theory is explained in detail. 



For a = 1, <5i is the energy associated with the discrete mean curvature operator discussed e.g. in [45[, while 
(?2 agrees with the quadratic form associated with and in fact the subdifFerential of <^2 is just the discrete 

Laplacian. If G is uniformly locally finite, then by Lemma 12.71 X and hence A2 are bounded linear operators. 
Accordingly, the initial value problem associated with 

Lp{t,v) + A2ip{t,v) ^0, felR, veV. 

is well-posed in ^^(V) and the exponential matrices e*^^, t e R, yield its solutions. Fiedler has showed that the 
largest eigenvalue of A2 is larger than the maximal node degree of G. Thus, dropping the assumption of uniform 
local finiteness A2 becomes an unbounded operator whose spectrum is not contained in any left half-plane, hence 
not the generator of a strongly continuous semigroup on ^^(V). Still, it can be proved by methods based on 
quadratic forms that — A2 is a generator. Our aim in the next section is to show that, just like in the special 
case of p = 2, most properties of the continuous p-Laplacian carry over to its discrete counterpart, despite the 
latter being non-local. 

3. Well-posedness results 

In view of the analogy between discrete and continuous calculus, the correct discrete pendant of the continuous 
p-Laplacian is the nonlinear operator on defined by 

(3.1) Ap/:=l(|I^/r2z^/), /GR^, 

in the case of an unweighted graph. More generally, we are going to consider the nonlinear operator given by 

(3.2) Cpf := (a|Z^/r-2(X^/)) , / £ R^. 

Definition 3.1. Let G be a graph. A family of graphs (G„)„gig is called growing if for all n,m G N with n < m 
G„ is an induced subgraph of Gm ■ It is said to exhaust G if for all n G h\ 

• G„ is an induced subgraph of G, i.e., G„ ~ (V„, E„,/i„,;/„) with /i„ = /ij^ and Vn ~ / o,nd 
•U„eNV„=V. 

The following is the main result of this section. 

Theorem 3.2. Let T > 0, fo G w^'^/iy) and f G L'^{0,T;ej{y)). Then the Cauchy problem 

, . ( Mt,^) = -Cp(bit,y)+fit), tG[0,T], vgV, 

^ ^' I 0(0, v) = /o(v), VGV, 

admits a unique solution cj) G iJ^O, 7"; ^d(V) H i°°(0, T; u;^'^'^(V)). 
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Furthermore, there is a growing family of graphs (G„)„giH that exhausts G and such that the sequence of 
solutions (0n)neiM to the Cauchy problen^ 

( 0„(t,v) = -/:^0„(i,v) + /(t), t e [o,r], vG v„, 

(HEp(")) <^ 0„(t,v) = 0, t e [0,r], vg7V(V„), 

[ </'t,(0,v) /o(v), V e V„, 

converges to 4>, weakly in -ff^(0, oo; ^^(V)) and weakly* in L°°(0, oo; t«^'2'^(V)). 

// moreover / = and fo G ^d(^) /'^'^ some q G [l,oo] (resp., /g G Cod(V)j, i/ien G C(IR+;£^(V)) (resp., 
(/> G C(1R+; Coc;(V) j and ('/>n)neiM converges to 4> also weakly* in i°°(IR_|_; i'^(V)) (resp., in i°°(lR+; cod(V)) j. 

Thus, we can approximate a dynamical system on the whole G by a sequence of dynamical systems on induced 
subgraphs G„. Before proving it, we will need a number of results. 

Lemma 3.3. Let p G (1, cxd). Then the proper functional S'p is convex and S'p + ^\\ ■ \\'^2 is coercive for any lj > . 
Furthermore, Sp is continuously (Frechet) differentiable as a functional on w^'2'^(V) and lower semicontinuous 
as a functional on ^^(V). 

Proof. Since for all p G (1, oo) 

(3.3) ^p = i|I.||^.(,)OX^ onwl'^/iW), 



the composition of a convex and a linear mapping, the functional S'p is convex - hence so is any perturbation by 
another convex functional, and in particular iS'p{-) + f || • ||^2 is convex and coercive for any uj > 0. 

In view of p.3p . and because l'^ is a bounded linear operator from ^^'^'^(V) to £q(E), in order to check 
continuous differentiability of S'p it suffices to observe that the functional || • ||^p on £p{E) is continuously differ- 
entiable for p G (l,oo). Thus, S'p is in particular lower semicontinuous as a functional on ^^'^'^(V), and lower 
semicontinuity as a functional on £^(V) follows from [s^, Lemma IV. 5. 2]. □ 

Remark 3.4. One can see that if additionally G is uniformly locally finite, then Sp is lower semicontinuous also 
for p = 1. However, we cannot invoke again fs^ . Lemma IV. 5. 2], as this result relies upon reflexivity of the 
effective domain. Nevertheless, by Lemma \2.T\ i^ is a bounded linear operator from ^'^(V) to ^^(E). Moreover, 
lower semicontinuity of \\ ■ ||£1(e) (viewed as a function from £ai^) to [0,+oo]) follows from Fatou's Lemma. In 
view of 

we also deduce lower semicontinuity of Si . From now on, we will always restrict to the case of p > \ . 
Lemma 13.31 yields that Sp has a subdifferential, which is at most single- valued by Lemma 



Lemma 3.5. Letp G (1, oo). Under the Assumptions^^ the subdifferential of Sp with respect to is given 

by 

^p/(v) i E - fHr^ifi^) - /W). V G V, 

^ ' wev 

w~v 

with maximal domain. 



Here and in the following N{\N) denote the neighborhood — in the graph theoretical sense - of the subgraph of G induced by 
some node subset W C V: i.e., JV(W) is the set of all v S V such that fi{v, w) ^ for at least one w G W. 
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Proof. Let / G ^^(V). If p G (l,oo), by Lemma [2.61 and because (op is Frechet differentiable whenever restricted 
to w^'^'^(V), its subdifFerential dS'pf C ^^(V) is either empty or a singleton. By virtue of 32, Lemma 2.8.8], for 
/ in the domain of Cp the only element of Cp is given by 

d^pf = g ^ ^;{f)q^^{g\cj,)^, V0G^2(v) 

^ ^a(e)|(I^/)(e)r2(I^/)(e)(I^0)(e)^5].g(v)0(v)d(v) Vcb e £j{y) 
eeE vev 

^ -i^l(a|X^/r2(X^/))(v)=.g(v), v £ V, 
in accordance with (j3.2p (recall that by assumption (i(v) > for all v G V). More explicitly, for all v G V 

p-2 / 



54./(v) = ^E^ve(a(e) 



wev 



X! ''we/(w) 



Vwev 



5] .ve («(e) |/(eo) - f{e,)r' {f{e,) - /(ei))) 

eeE 



^^-^ wev 

This completes the proof. □ 

Proposition 3.6. Let p G (l,oo). The set of eigenvectors of Cp for the eigenvalue is a subspace of £^{\/) 
spanned by the characteristic vectors of the connected component of G with finite surface, in the sense of Defini- 
tion\£E 



Proof. If f G ^d(^) eigenfunction of Cp with eigenvalue 0, then f = 0. By Remark l2.8l this implies / = 
if and only if each connected component of G has infinite surface. □ 

By Lemmas 13.31 and 13.51 a direct application of [s^ Prop. IV. 5. 2] yields the following. 

Proposition 3.7. Let p G (l,oo) and let a,d satisfy Assumptions \2.4\ Let T > 0, /o G and f G 

L2(o,T;£^(V)). Then the Cauchy problem (HEp) admits a unique solution (j) e H\0,T;£l{Vj)nL°°{0,T;wl-P/{Vj). 

If f then the solution is given by a strongly continuous semigroup of nonlinear contractions and cj){t, ■) G 
D{Cp) for all t > 0. 

In order to prove relevant properties of the p-heat equation one often wants to prove that certain relevant 
subsets of the state space are left invariant by {e~*'-'p)t>o. This can be done by a nonlinear generalization of the 
Beurling-Deny criteria due to Barthelemy In the continuous case most proofs are made easy by the locality 
of the operator. In the discrete case locality fails to hold: The key point in our proofs will then be that the 
considered orthogonal projections map vectors in ^^(V) into vectors with smaller oscillation. 

We will need the following apparent property of the ^^'-norm. 

Lemma 3.8. Define a function 

fk^p : [0,oo) 3 a^\k + a\P + \k ~ a\P G [0, oo). 

Then fk.p is strongly monotone increasing for all k Cz R and p G [1, oo). 

The following is a discrete analogue of various results proved or summarized in [ill . § 4.1]. In the special case 
of p = 2, these results have been obtained in [g^] already. 
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Proposition 3.9. Let p G (l,oo ). T he strongly continuous nonlinear semigroup (e *^'')t>o generated by —Cp is 
(sub)Markovian, in the sense of \1 i . Def. 2.2.], i.e., it is order preserving, 

f <9 e-^^^f < e-^^^g Vi > 0, 

and II • -contractive, 

||e-*^P/„e-*^''g||,~ < ||/-.g||,~ > 0. 

Moreover, it is positivity preserving, 

/ > ^ e-^^^-f > Vi > 0. 
Proof. Because S'p is nonnegative, the eondition 

(3.4) A g) + Sp{f V g) < Sp{f ) + Sp{g) for all /, g e w]:f{y), 

is sufficient (and necessary) for the semigroup to be order preserving, cf. 0, Cor. 2.2]. In particular, the semigroup 
is order preserving if for all e G E and all /, g G w^'^' (V) 

|X^(/ A g){eW + |X^(/ V g){eW < |I^/(e)r + l^'^gie)]"- 
This can be proved taking e G E and /, g G ^^'^'^(V) and dividing the four possible cases 

• /(eo) < .9(eo) and /(ei) < g{ei), 

• .g(eo) < ./(eo) and ^(ei) < /(si), 

• .g(eo) < /(eo) and /(si) < .g(ei), 

• /(eo) < .9(eo) and g{ei) < /(ei). 

The assertion clearly holds in the first two cases. In order to check its validity in the third and fourth case, we 
prove that for all x,y,w, z G R with x > y and w > z 

\y - z\P + \x- w\P <\x- z\P + \y- w\p, 

i.e., 

\k + a\P + \k- a\P > \k + I3\p + \k - f3\P, 

where 



X + y — w — z X — y + w — z 
k := z , a := , p := 



-X + y + w — z 



2 2 2 

In fact, ii X > y and w > z, then 

y — x + w — z<y — X — w + z<x^y^w + z, 

i.e., |a| > 1/3 1 and the assertion follows from Lemma \3M The other cases can be treated likewise. 

Let us now prove that the semigroup is || • H^t^ -contractive. First of all, because S'p is homogeneous (of degree 
p)j by [ll|, Cor. 3.7] we have to prove that 

'.9 + (/-,9 + l)+ , g-if~g~l)- 



2 2 

/-(/-.9+l)+ , /+(/-.9-l)- 



2 2 



< Sp{!) + Sp{g) for aU /,5 G w^fi^^M). 



As above, it suffices to take one e G E and /, g G w^p^ (V) and to check that 

(e)+2: (e) 



'/-(/-.9+1) 



±(e)+I^^ + (^'g'^^- 



2 ^ ' 2 

< |X^/(e)r + |I^g(e)r for all /,g G 



■(e) 
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or equivalently 



z + [w — z + 1)+ y — (x — y — 1)_ z — — z — 1)_ 



{x-y + l)+ 



{w — z + X + (x — y — 1)_ w + {w — z 



2 2 2 

< \x - w\P + \y - z\P 

for all a;, y,w, z G R. Proving this inequality in the nine possible cases 
y\ < 1 and |w 



1) " 
-^(ei) 



• \x — y\ S and |w — z| < 1, 

• ^ 2/1 !i 1 E^nd w — z > 1, 



y\ < 1 and w — z < — 1, 



• 


X 


- y 


< — 1 and \w 


-z\<l, 


• 


X 


- y 


< — 1 and w 


- z > 1, 


• 


X 


- y 


< — 1 and w 


- z < -1, 


• 


X 


- y 


> 1 and \w — 


z\ < 1, 


• 


X 


- y 


> 1 and w — 


z > 1, 


• 


X 


" y 


> 1 and w ~ 


z < -1, 



is tedious but not difficult. E.g., in the second case (the first being trivial) one has to check that 



w ■ 



w + z — 1 



<\x- w\P + \y - z\f- 



Also in this case, this condition can be re-written as 



where again 



and 



-o\P > |fc + 7|P + |A:-7|f, 

z —X + y + w — 



7 



y 



1 



Under the assumption that (|a; — j/| < 1 and) w — z > 1, one has 

X — y — w^z<x — y^\ < y — x 



and hence I7I < \a\, whence the assertion follows, by Lemma [3.81 

Finally, by 0, Cor. 15.5] the semigroup is positivity preserving if and only if 

<^p{n < Sp{f) for all / e 
This is actually satisfied, because for all e G E dividing the cases 

• /(eo),/(ei)>0, 

• ./(eo),/(ei) <0, 

• ./(eo) < < /(ei), 

• ./(eo) > > /(ei). 
one sees that 

|/+(eo)-/+(ei)|<|/(eo)-/(ei)|. 

This yields the claim. 



□ 



Remark 3.10. In fact we have just showed more than we have stated: our proof yields that the same assertions 
hold for any restriction of Cp, in particular for any restriction to induced subgraphs. 
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Remark 3.11. Let S he a coercive, convex, lower semicontinuous functional on a Hilbert space H. By 
Prop. 4.5], if a closed convex subset of H is left invariant under {e~^^'^)t>o for all t > 0, then the same holds 
under all resolvent operators J\ := (Id + \dS)~^ for all A G (0, 00). 

1) In this way one easily sees that e~*^p is not homogeneous (of any degree) for any p ^ (l,oo) and anyt>Q, 
although Cp is homogeneous of degree p. However, e~*^pO because S'p{0) = 0. 

(2) In particular. Proposition \3.9\ yields the following comparison result: 

Let A > 0. If f > 0, then the solution ip of the elliptic problem 

(3.5) A(^(v) + Cpip{\/) = /(v), V G V, 

satisfies ip >0. Moreover, if ipi,ip2 are the solutions of the elliptic equation with inhomogeneous data /i,/2 and 
if fi ^ /2; then also > 1^2- This result is comparable to results that are known to hold on finite unweighted 
graphs, cf. fi3 /. 

Lemma 3.12. Let p G (1,00). The strongly continuous nonlinear semigroup (e~*^i')t>o extrapolates to a family 
of nonlinear semigroups on ^^(V) for all q G [l,C)o] (strongly continuous if q^ 00) as well as on 

Cod(V) := {/ : V ^ R : lim /(v„)rf(v„) = 0}. 

n— >cxD 

Formally speaking, this result yields well-posedness of a dynamical system only if one is able to determine the 
domain of Cp. In the linear case (p = 2) this has been done in (27j . 

Proof. By Proposition l3.9l contractivity of (e~*'^'')t>o with respect to the norm of yields that the semigroup 

on ^^(V) extends to a contractive semigroup on the closure of ^^(V) in the ^^(V)-norm, i.e., in co(i(V). By duality 
we obtain a contractive semigroup on i\(y), and finally again by duality a contractive semigroup on l^iy). Now, 
the assertion follows applying the non-linear Ricsz-Thorin-type interpolation theorem due to Browdcr, see [ill . 
Thm. 3.6] for a version tailored for our setting. □ 

Proof of Theorem VJ. 2\ As already emphasized in Remark I3.10[ the proof of Proposition 13.91 shows the [| • ||£<=o- 
contractivity of the semigroup generated by the subdiffcrential of each restriction of Sp. Consider the sequence 
of finite dimensional subspaces of ^^(V) of the form 

^2(V„) = {/ :G f^(V) : /(v„0 = Vm > n} = RV", 

where V„ := {vi, . . . , v„}. Accordingly, the restriction of to ^^(V„) is the functional given by 



<"'/W^=i E «(v,w)|/(v)-/(w)r + i «(v,w)|/(v)|P, vgV. 

vev„ 

veAr(V„) 



^ v,wev„ ^ vev„ 



A direct computation yields that — £p"'' := dSp^^ is the restriction of —Cp to the space of function on the 
subgraph induced by V„ U N{\/n), with Dirichlet conditions on 7V(V„). Hence, the associated Cauchy problem 
is (HEp*^"'). Fitting to the present setting the usual arguments that yield on convergence of the Galcrkin 
scheme (e.g., following the proof of ^ Thm. 6.1, Part 3]), one sees that the sequence of solutions of (HEp(")) 
converges (up to subsequences) to the solution found in Proposition [221- weakly in i/i(0, T; £^(V)) and weakly* 
in i~(0,T;«;^;^'2(V)). 

Let now / = 0. Then for all n G N the solution to each (HEp'")) is given by a strongly continuous semigroup 
of nonlinear contractions, and the same holds for the the solution to (HEp). Because the restricted energy 
functional S'p^^ is a proper, convex and lower semicontinuous functional, (minus) its subdifferential generates a 
semigroup on ^^(V„); which is || • ||ft^ -contractive by ProDOsition l3.9l By the same interpolation argument used 
in the proof of Lemma r3.121 this semigroup extrapolates to cod(V) as well as to ^^(V) for all 9 G [l,oo]. In 
particular, for all /o G ^^(V) 



" /o||i=(0,oo;£«(V„)) < ll/0|v„ ll<?;j(V„) < ll/oll<?9(V)i 
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due to the fact that (S'p(O) = and hence e p = 0, and because of the contractivity of the semigroup. Because 
for all q G (1, cxd] L°°(0, oo; ^^(V)) = L^{0, oo; £^ (V))' is the dual of a separable Banach space, by the theorem of 

Banach-Alaoglu for all < > the bounded sequence {e'^^p '/o)neN C L°°(0, oo; £^(V)) has a weak*-convergent 
subsequence. □ 

Recall that a semigroup on a Banach space X is called irreducible if the only ideals of X that arc left invariant 
under it are the trivial ones. 

Proposition 3.13. Let p G (1,cxd). Then the graph G is connected if and only if the semigroup (e~*^J')t>o is 
irreducible. 

In the special case p = 2, by [s^ Thm. 2.9 and Def. 2.8] connectedness hence implies that if / > but / ^ 0, 
then e~*^^f > 0: this property of the discrete Laplacian is already known, see e.g. 27|, Cor. 2.9]. 
The proof of sufficiency is due to Rene Propper (Ulm). 

Proof All ideals of £^(V) are of the form ^^(Vo) = {/ G £^(V) : supp/ C Vq} for some subset Vq C V, and the 
associated orthogonal projections are given by the restriction operators Pvo ■= Ivo'- 
Now, assume ^^(Vq) to be invariant under (e~*'^p), or equivalently that 

(3.6) '^piPvof) < for all / G £^(V). 

We have to show that Vq is a trivial subset of V, i.e., Vq = V or Vq = 0. 

In fact, if Vq 7^ V 7^ Vq^ there are two adjacent nodes vq G Vq and vi G Vg'. Set 

V:=V\{vo,vi}, 

so that E is partitioned into (vg, Vi), E°, E^, E, where for i G {0, 1} E' consist of those edges other than (vo,Vi) 
one of whose endpoints is (regardless of their orientation), and E := E \ ((vq, Vi) U E" U E^). In other words, 

pS'pig) = a(vo,Vi)|.g(vo) -.g(vi)|^' + ^ a(vo,w)|g(vo) - g(w)|'' 

W~Vo 



+ J2 a(vi,w)|.9(vi)-.g(w)r + ^a(e)|.g(eo)-g(ei)r for all g G ^^/(V). 
Let now / G ^^(V) be defined by 



v#vo 




if V = vo, 
if V = vi, 
0, otherwise, 



for some a; > to be determined later. Accordingly, 

^VoMvj • I otherwise. 

Therefore, 

pS'pif) = a(vo,Vi)|a; - 1^ + ^ a(vo,w)|a-|P + ^ a(vi,w). 



(observe that the sums on the RHS are finite, because they are less then deg(vo) and deg(vi), respectively 
recall that G is assumed to be locally finite) whilst 

pSp{PyJ) = a(vo,vi)|x|f + a(vo,w)|x|P. 

W~Vo 
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Accordingly, 

pMf)^PMPvof) = a(vo,Vi) (|x - l\P - \xn + J2 «(vi,w) < 0, 

v#vo 

choosing x large enough, thereby contradicting p.6p . This implies that indeed Vg = or Vq = G. 

Conversely, each subspace of ^^(V) consisting of functions over a connected component is apparently left 
invariant under (e~*^f )t>o- □ 

4. Symmetries 

Graphs bearing some symmetry are well-studied objects of graph theory. In this section we are going to show 
how non-trivial symmetries of (HEp) arise if the underlying graph enjoys special symmetry properties. Because 
S'p and Cp do not depend on the orientation of G, we can fortunately neglect it throughout this section. 

To begin with, we need a simple lemma. 

Lemma 4.1. Let Hi,H2 be two Hilbert spaces and Ei : Hi — >■ [0,oo] and E2 : H2 — !• [0,cxd] two proper, lower 
semicontinuous, convex functions. Let T, be a bounded linear operator from Hi to H2 ■ Then E intertwines with 
the semigroup generated by —dEi and — 9i?2, i.e., 

g-taiJi ^ ^ ^e-tdE2 allt>0 

if and only if 

El {Lf + j:*Rg) + E2 {^Lf + g- Rg) < Ei (/) + E2 (g) for all f e Hi, g e H2, 

where 

L := (/h, and i? := (7^, + SS*)-^ 

Proof. One checks directly that S intertwines with the semigroups if and only if the graph of S, i.e., the closed 
subspace 

Graph(S) := | (^.^^ e Hi x H^ 
is invariant under the matrix semigroup 

/„-tdEi n 

-tdE I e U 



e~*^^= ; ' ^ - °' 

where E :— Ei ® E2. A classical formula due to von Neumann yields that the orthogonal projection of Hi x H2 
onto Graph(E) is given by 

■ L 



^Graph(E) - yj.^ _ ^ 

cf. [H Thm. 23]. The assertion follows from [3, Theo. 1.1]. □ 

Definition 4.2. A permutation O on\/ is called a node automorphism of G ~ (V, E,/x, z^) if for all v,w £ V 

• v{Oy) = J^(v) and 

• the entries of A introduced in (|2.1|) satisfy aov Ow = Q^vw (i.e., {0\/,0\n) G E or {0\N,Oy) G E if and 
only if {\/,w) € E or (w,v) G E, and in this case ^{Ov,Ovj) = /i(v,w)j. 

We denote by Aut(G) the group of all node automorphisms of\/. 

In the unweighted graph, the above definition reduces to the usual one: node automorphisms arc permutations 
on V that preserve the adjacency relation. We emphasize that this definition docs not depend on the orientation 
of G. 

Observe that each node automorphism O induces a permutation Ol on E by 

OLe~{0\/,Ow) ife=(v,w). 
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Thus, each edge permutation O defines a mapping on and a mapping on defined as the Nemitskii operators 
associated with O and Ol, respectively, i.e.. 



and 



u M> u{Ol-)- 



With a shght abuse of notation, we wiU not distinguish between the node/edge permutations and their associated 
Nemitskii operators. 



Theorem 4.3. Let p € (l,oo). If O E Aut(G) and a = fi as well as d=v (i.e., Cp = Ap 



), then 



-,-tA 



for all t > and fc G N. 



If G is unweighted, this resuh is clear for p ~ 2, since then A = D — A ioi a. diagonal matrix D. 

Proof. Clearly, it suffices to prove the claimed commutation relation for k = 1. This can be checked owing to 
Lemma [41] observing that setting Y, := O one has 



Indeed, 



f + 0*g 



and 



(Z^/)(e) + (I^O*.9)(e) 



(I^g)(e) + (J^O/)(e) 



By a change of variables and using the fact that a is constant along induced orbits we obtain 



f + 0*g 



OieeE 
1 1 



Of + 9 



(I^O/)(e) + (T^g)(e) 



eeE 



(T^O/)(e) + (J^g)(e) 



p2P- 



jY.a{e)\{I^Of){e) + {I^g){e) 



eeE 



< - fE«(^) (|(^^0/)(e)r + |(I^.9)(e)f )) 

^ VeeE / 

= i('^a(e)(|(I^/)(e)r + |(I^5)(e)r 

^ VeeE 



<^p(/)+4>(5), 



as we wanted to prove. 



□ 



Remark 4.4. (1) Let G have finite surface. Averaging a function u over all nodes (i.e., shorting all nodes, in 
the point of view of Remark \2.10\) one obtains a system which is trivially left invariant under the evolution of 
the p-heat equation, since the associated projection P satisfies of course (S'p(l) = 0. This is independent on the 
automorphism group of Q, see Remark \4. 10\ (2) below. 

(2) On the other hand, just shorting two arbitrary nodes is not sufficient to obtain an invariant subsystem: 
this can be easily seen by taking a path of length 3 with /i = 1 and considering the projection P of (V) = R^ 
onto the space 

{u : {vi, V2, V3, V4} R : u(v2) = ""(vs)}. 
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Then, for m(v„) := n one has 

which by f^, Theo. 1.1] shows that said subspace is not invariant under e~*'^p. 

Clearly, each subgroup of Aut(G) defines a partition of V into equivalence classes with respect to their orbits. 
We denote by [v] such orbits and by \[y]\d their lengths with respect to d, i.e., 

MU — Yl 

we [v] 

Theorem 4.5. Let p G (1,cxd). Let T be a subgroup o/ Aut(G) with orbits of finite length with respect to v. Let 
a = fi as well as d = v (i.e., Cp — Ap as in (|2.2p ) and consider the orbitwise averaging operator P defined by 

PfM ■■^tAtT. /(w)Kw), / e iUy), V e V. 

Then the following assertions hold. 

(1) The range of P is left invariant under {e~*'^'')t>o- 

(2) The null space of P is left invariant under {e^*'^^)t>o if P = 2. 

Example 4.6. 1) A typical and relevant case is that of an infinite radial trees T, i.e., infinite rooted trees such 
that any two nodes with the same distance n from the root have the same number 7„ of children - for the sake of 
simplicity, say, T is binary, i.e., 7„ = 2 for all n £H. In this case, the automorphism group ofT is abstractly 
isomorphic to a semidirect product of permutation groups and any orbitwise projection acts by simply averaging 
-functions supported in an infinite binary subtree Tq o/T over all nodes at the same distance from the root of 
Tq. Then, Theorem \4.5\ shows that radial initial data give rise to radial solutions. This is already known in the 
linear case p = 2, both in the case of discrete and metric graphs f fsdl . Udll). The case of graphs that are not trees 
has been studied in fl^ l. 

2) The range of P with respect to the subgroup T is isomorphic to a new "quotient graph" G/T obtained 
identifying all the nodes belonging to the same orbit; in the case of radial trees, if T is the full automorphism 
group, then G/T is a semi-infinite path, i.e., N. 

Instead of proving the above theorem directly, we will deduce it as a corollary of a more general result. 

Definition 4.7. Let I be a (possibly infinite) set. A partition 

v = Uv. 

of the node set \/ of G is said to be almost equitable with cells {\/i)ii=i if for all i,j £ I , i ^ j, 

(4.1) there are numbers Cij > s.t. //(v, w) = Cyi'(v) for all v g V^. 

w6Vj 

We write 

|V,|^ := ^ i^(v), iel. 

Remark 4.8. (1) Clearly, each node partition of a graph canonically induces an edge partition: Simply take 
edge cells E^j consisting of all edges with initial endpoint in \/i and terminal one in Vj ( orientation does matter! 
- that is, Eij 7^ Eji)- Moreover, we let [e] := if e £ Eij and write 

(4.2) |Ey|a := '^(^)' ^ P 
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(2) In particular, if is an almost equitable partition, then = Cij\\/\u for all i,j G / s.t. i ^ j. 

(3) If the condition (|4.ip holds for all i, j (and not only for i ^ j), then the partition is called equitable. 
Our definition is a generalization of the classical one for unweighted graphs (fi = 1, v = \), see e.g. f2dl. Isdll for 
equitable partitions and 0, § 2.3] for almost equitable ones. In the unweighted case, existence of an equitable 
node partition amounts to saying that each node in is adjacent to exactly Cij nodes in Vj . 

Then, the quotient graph is the directed, weighted multigraph with node set of cardinality \I\ (the i-th node 
\Ni corresponding to the cell V^J such that (wi,Wj) is an edge (and if so, with weight Cij) if and only if Cij ^ 0. 
A simple but useful result in the linear, finite, unweighted case is that the spectrum the discrete Laplacian of G 
contains that of the discrete Laplacian of its quotient graph, see e.g. Jsdi . Thm. 2.3]. In certain special cases the 
spectra even agree ( not counting multiplicity, of course ), cf. Jldi . Thm. 7. 8]. 

Theorem 4.9. Let p G (l,oo). Let G have an equitable partition associated with cells s.t. \\/i\d < oo and 

\^ij\a < OO for all i,j £ /. Let a = fi as well as d = v (i.e., Cp = Ap as in (|2.2p ). Assume that 

there exist aij > s.t. a(e) = atj for all e G Eij, and all i,jEl, i =/= j. 

Consider the cellwise averaging operator P defined by 

PfM E /W^W. / e ^d(V), V e V,. 

Then the following assertions hold. 

(1) The range of P is left invariant under (e~*^p)t>o. 

(2) The null space of P is left invariant under {e~*'^'')t>o if P = 2. 



This can be seen as a nonhnear counterpart of a well-known result of algrebraic graph theory, see e.g. [21 
Thm. 9.3.3], or a discrete counterpart of the settings in [sl. l4o| (where different terminologies are used). 

Proof. By [H, Theo. 1.1], the range of P is invariant under (e~*'^'')t>o if and only if 

(4.3) \\:i''P.m. < Wl'^m. for all / E w'-fiy), 
whereas its null space is invariant if and only if 

(4.4) ||X^(/ - PfW,. < III^/II^P for all / e w^^^'iv). 
First of all, define the averaging operator 

Pu{e) "(0«(f), « G ^'(E), e S E s.t. e G Ey, 

and observe that for all / G ^^(V), all e G E, and all i,jEl 

• I^P/(e) = if e G E;,, since P/(v) = P/(w) for aU v, w G V;, 
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while if e G with i ^ j, and hence bq G Vi and ei G Vj, then 



|P/(eo)-P/(ei)| 

1- ^ /(v)d(v) - ^ E 



IV 



veVi 



wev, 



(*) 



|v,u 

alE 



P(V) 



-ij I a 



f6Ei 



where (*) follows from the definition of equitable partition and Remark 14.81 f2). 

(2) This shows that, for p = 2, we are actually looking for properties of the variance and the expected value 
of the random variable X := f: hence it suffices to apply the formula Var{X) + (£'[X])^ = E\X\'^ to in order 
to check that the claimed inequalities hold for p = 2. 

(1) In order to prove (|4.3p for general p > 1, we proceed as follows. One has 



eeE 



iG/ eGE;, 



E a. E \P^^fi^)\ 



< 



i,jel eeEij 



E '^'^ E 
E E 



I fGE.,- 



fGE, 



where the last step follows from Jensen's inequality. In view of (|4.2p . wc obtain 



^a(e)|Z^P/(e)r < E E l^"^ ^ E 



eeE 



i,je/ feE. 



eeE 



This shows that (j4.3|) holds and completes the proof. 



□ 



Remark 4.10. (1) The second assertion in Theorem \4-.5\ is optimal, in the sense that there exist graphs on which 
the null space of P is not left invariant under e^^^^ for any p ^ 2. An example for which (|4.4p fails to hold 
for any p ^ 2 is given as follows: Take G to be a path of length 3, with weights fi = 1 and v = \. Consider the 
trivial equitable partition := {v} for all nodes) and the function 

1 



/= ( l,-,0,x) G 
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for X ^ R to be fitted later. Then 



and choosing values of x slightly larger than 1 for p > 2, and slightly smaller than 1 for p < 2 yields the 
sought-after counterexample, by elementary calculus arguments. 

(2) We have already observed in Remark \4.4\ (V ihat averaging over all nodes of a graph with finite sur- 
face yields a projection whose range is invariant under {e~*^'')t>o for any p. This is in fact a special case of 
Theorem \4-.9\ since each graph has the (trivial) almost equitable partition given by Vi = V. 

(3) Theorem \4.9\ is strictly more general than Theorem \4-.5\ since the orbit partition w.r.t. any subgroup T of 
Aut(G) yields an equitable partition of\J, but the converse is generally false (cf. f2dll). 

(4) The (non- equitable) partitioning of\/ into Fix(T) := {v G V : Ov — v VO S F} and its complement is in 
general not respected by {e~^^'')t>o for any p, as the simple following example shows: Take Gm as an unweighted 
m-star with a path of length 2 attached to the center, and consider the function f defined by f{v) = on each 
node V of the star (including the center) and f{v) — 1 on both the remaining nodes. 

5. Appendix: Further discrete operators 
In this section wc comment on several different possible extensions of the theory introduced above. 

5.1. Generalized Laplacians. Fiedler, Colin de Verdiere and others have introduced and studied several vari- 
ations on the classic discrete Laplacian. The following one, which goes back to jl^l, is the most usual one: Given 
an simple, finite, connected graph with node set V, any |V| x |V|-matrix whose off-diagonal v-w-entry is 

• = if and only there is no edge between v and w and 

• < otherwise 

is called a generalized Laplacian. Besides the discrete Laplacian, also —A (where A is the adjacency matrix) and 
Chung's normalized Laplacian 

(E § 1-4]) 

V := D-^AD-i 

(where D is the diagonal matrix whose i-z-entry is the degree of v^) are clearly generalized Laplacians. Unlike 
the adjacency matrix, Chung's V is also positive definite, since 

P = {D-^I){D~^lf . 

The signless Laplacian 

Q:=D + A 



has been introduced in |14| . Its study has gained much momentum in the last decade, also owing to thorough 
investigations by Cvetkvovic and others, cf. [l^j. (Clearly, — Q is a generalized Laplacian, too). One possible 
reason for this popularity is the richness of interesting graph theoretical properties of Q (e.g., is always an 
eigenvalue whose multiplicity is the number of bipartite components); another one is its nice variational structure. 
In particular, one easily sees that 

Q JJ^, 

where is the incidence matrix of the undirected graph underlying G, i.e., J" :~ 1+ +T^. One may further 
generalize this class of Laplacians by taking a G £^ (V) and letting 

( J,)ve X+ + a{v)I-^ , V e V, e e E. 

We may then introduce 

which is still a positive definite matrix and also a generalized Laplacian (each of whose non- vanishing off-diagonal 
entries is one of the CTv); along with its normalized version 
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(This generalization process could continue by allowing for more general matrices a, but this would necessarily 
destroy locality of 2°^). It is then natural to consider a signless p-Laplacian by 

Qpf:^J{\J^.fr'sign{J^f)), 

or, more generally, 

Q^/:= J.(|jJ/r^sign(jJ/)). 
Likewise, one can introduce the normalized operator 

P;f D-iM\j;fD-ifr'sign{jjD-'^f)). 

We are aware of only a few previous investigations on the signless p-Laplacian in the literature, including [3l[23j. 
where it is proved that the sets of eigenvalues of Qp and Ap agree, provided that G is bipartite. 

Remark 5.1. In the linear case p ^ 2, the parabolic theory of the signless Laplacian on uniformly locally finite 
graphs is not overly interesting. This is due to the fact that Q + A = 2D (we are omitting the index 2 o/ A2, Q2)- 
All these three operators are bounded. Although D and A (and therefore D and A, do not commute, so that 
in general 

e*2 ^ e-^'^e^*^, t £ R, 
Lie 's product formula still holds and yields 



e*2= lim (e-i^e^i^Y, te 

n— ^00 V / 



Alternatively, we we can recover another formula for {e^^)t>o based on the Dyson-Phillips formula. Since 
D is a positive diagonal matrix, knowing qualitative information on (e^*^)tgR promptly yields quite complete 
information on (e*2)tgR, and vice versa. For example, since (e~*^)t>o is a positive semigroup, so is (e*®)t>o. 

Things are different in the general nonlinear case, since it is not clear in which sense Qp can be seen as a 
perturbation of Ap ( unless p is an even natural number, in which case Ap and Qp can be compared using the 
binomial formula). 



The following can be proved like in the case of Ap. 
Proposition 5.2. Let p G (1, 00) and consider the functional ^p : ^^(V) — > [0, 00] defined by 



eGE 

Then the following assertions hold. 
(1) The effective domain 



z 



a,d 



(V) :={/e^^(V): jJ/e^g(E)}. 



of ^p is a Banach space with respect to the norm defined by 



\\f\\% + \\jnre- 



For all p G [l,oo], -Zq'^'^(V) is continuously and densely embedded into •^^(V). If moreover p € [1,cxd), then 
^a'd '^O^) separable. If p G (l,oo), then ^^'^'^(V) is uniformly convex (and hence reflexive). 

(2) The functional ^p is proper and convex. It is continuously (Frechet) differentiable as a functional on 
z^'^'^(V), while it is lower semicontinuous as a functional on 

(3) The subdifferential Qp of ^p generates a strongly continuous semigroup of nonlinear contractions on £^(V). 
(4-) If <J = 1, then the set of eigenvectors of Qp for the eigenvalue form a subspace of £^{\/) whose dimension 

is the number of bipartite components of G with finite surface. 

Observe that unless u = ±1, ^p and hence its subdifferential Qp do depend on the orientation of G. 
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Proof. Because 



and because is a bounded linear operator from 2^ d'^(V) to ^^(E), all the proofs are analogous to those for 
the corresponding assertions about and Ap. The only minor change is needed in the proof of (4), where the 
dimension formula for the null space of has to be replaced by a dimension formula for the null space of J^^ , 
cf. [13). □ 

5.2. The p(e)-Laplacian. Wc could also consider the equivalent of the so-called p(x)-Laplacian: in our discrete 
setting the associated energy functional is given by 

eeE ^ 

under suitable assumptions onp G R^. In this case it would not in general be the subdifferential of d'p, but rather 
that of its convex, lower semicontinuous relaxation that generates a semigroup. 

5.3. Discrete Schrodinger operators. Another possible direction of generalization consists of considering an 
energy functional with additional terms defined by 

i ^ a(e)|/(eo) - /(ei)!^ + ^ ^ 
^ eeE vev 

for some b asymptotically comparable with d. Due to their interpretation in (linear) potential theory, the terms in 
the second sum are sometimes referred to as killing terms. The subdifferential of S'p is in fact the discrete analogue 



of a Schrodinger operator with scalar potential b. We refer to 27[ and references therein for a comprehensive 
theory in the linear case. 

5.4. Discrete operators with boundary conditions at infinity. One can introduce the space i&q'^'^(V) 
defined as the closure of coo(V) with respect to the norm of w^'^'^(V). Whenever G is finite, u'|^'2'^(V) always agrees 
with iL'^'2'^(V) (in fact, both agree with £^(V)), but for infinite graphs things are less obvious. A characterization 
of the quotient space 

is a classic topic of potential theory and in special cases it leads to the introduction of the so-called Martin 
boundary (see e.g. [35| for its connection to the notion of space of ends of an infinite graph), but in the general 
case we are only aware of a subtle discussion of parabolic and graph theoretical properties implying non-triviality 
of u'^'^'^(G)/'u;^'^'^(G) performed in [ij, §§ 4-5]. While the potential theory of discrete boundary value problems 
seems to be confined to the linear case therein, one may introduce also in our case the quotient Banach space 

If Tr denotes the canonical surjection of i«^'^'^(G) onto b^ d(9G), one may consider the general functional 

:=^||I^/||?. + i||Tr/||^,^(,,,. 

Then, it is easy to check that S'm,q is proper, lower semicontinuous and convex and the associated subdifferential 
is a discrete p-Laplacian with discrete Robin-type boundary conditions. 
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5.5. The porous medium equation. We briefly discuss the interplay between porous medium equation and 
p-heat equation on graphs. For the sake of simplicity, we only consider the unweighted case, i.e., /i = a = 1 
and v = d = 1. The continuous porous medium equation can be studied in the context of the theory of 
subdifferentials, see e.g. [38|, Examples III.6.C and IV. 6. B]; but it is known that in the continuous, 1- dimensional 
case if tJj solves the porous medium equation 

then there is (j) such that tJj is its pressure (i.e., (j)x = ip) and it satisfies the p-heat equation 

4> = i\4'xf~'^4'x)x 



for p : / = TT + 1. While the proof in is based on the theory of exact differential forms (see e.g. [43|, § 3.4.3]), and 
it is not clear whether it has a pendant in our context, this argument roughly suggests that the correct space to 
discuss the porous medium equation is £^(E), which we can look at as the node space of the line graph of G 
and consider as a pressure space. We can prove a similar result in a special class of oriented bipartite graphs. 

Observe that bipartition induces a natural orientation: if V = Vi U V2, then assume all edges to have initial 
endpoint in Vi (and hence terminal cndpoint in V2). Wc recall that a bipartite graph (V = Vi U V2) is called 
semiregular if all nodes in have same degree di, i = 1,2. 

Proposition 5.3. Let G be uniformly locally finite and p > 2. Let T > and f G L^(0, T; £^(V)). If <j> ^ 
H\0,T;e(\/)) satisfies 

^{t, v) = -I(|I^0|f-2l^0)(i, v) + /(t), t > 0, V G V, 
then i) := (j) G H^{Q,T;P{V)) satisfies 

V'(<, e) = -I^XdV-r- V)(i, e) + X^/(i), i > 0, e G E, 
for p ~ TT + I. Furthermore, ip satisfies 

(5.1) e) = Ag, iW-'m, e) - (r + s)(|V^r- V)(t, e) + f{t), i > 0, e G E, 

if additionally G is (r^ s)- semiregular bipartite and is oriented accordingly. 

Observe that (|5.ip is a forward porous medium- type equation (with potential) on the line graph Gl of G. We 
have denoted by Ag^ the (linear) discrete Laplacian on G^. We stress that the latter result is not independent 
of orientation of G. 

Proof. It is a direct consequence of Lemma |2 . 71 that ip — i^4> if G is uniformly locally finite. 
Let now G be bipartite and oriented accordingly. The second assertion follows recalling that 

I^X = (r + s)Id - Ag^ 

whenever G is (r, s)-scmircgular, see e.g. [s^, Proof of Thm. 3.9]. □ 
Consequently, by Theorem 13 . 71 and Corollarv l2.9l the following holds, since an infinite graph has infinite surface 

Corollary 5.4. Let T > 0. If in particular G is a forest each of whose connected components is an infinite, 
uniformly locally finite tree, then for all g G L^(0, T; £^(E)) and all ipo G ^^(E) there exists a unique solution 
iP G iJi(0,r;£2(E)) to 



(PME) 



ij{t, v) = Z^I(|^|— ii^)(i, v) + g{t), t > 0, V G V, 

1^(0, v) = V'o(v), vgV. 
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